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Non-Abelian Anyons exist in certain spin models and may exist in Quantuam Hall systems at certain filling 
fractions. In this work we studied the ground state of dynamical SU (2) level-re Chern Simons non-Abelian 
anyons at finite density and no external magnetic field. We find that in the large-re limit the topological interac- 
tion induces a pairing instability and the ground state is a superconductor with d + id gap symmetry. We also 
develop a picture of pairing for the special value re — 2 and argue that the ground state is a superfluid of pairs 
for all values of re. 



Introduction. In two spatial dimensions, indistinguishable 
hardcore quantum particles are not restricted to being either 
fermions or bosons and can be anyons [llf]. The simplest 
anyons have abelian statistics, for which the wavefunction 
acquires a phase 6 different from zero or tt when the posi- 
tions of two particles are exchanged in a counter-clockwise 
manner, thereby implementing a form of topological interac- 
tion. Quasi-particles in many Quantum Hall states, such as 
v = 1/3, are expected to have such abelian fractional statis- 
tics 12 [3|. A finite density of such fractional quasiparticles 
can form its own fractional quantum Hall state, thereby lead- 
ing to states such as v — 2/5 (which also has anyonic quasi- 
particle excitations above its ground state) 10,3]. Anyons with 
abelian fractional statistics also appear in microscopic models 
of frustrated magnets in zero field , and one interesting 
question to ask is what is the ground state of a system with 
a finite density of anyons in zero field. It was found that the 
ground state of abelian anyons which do not interact with each 
other except through their topological interactions is a super- 
fluid IlZI |8|] . The physical picture is that the motion of a single 
anyon is frustrated since its phase gets scrambled due to the 
topological interactions with all the other anyons. Therefore, 
its kinetic energy is highly frustrated. A collection of anyons 
might, for some values of 9, behave like bosons and therefore 
may condense to form a superfluid. An argument based on 
a mean-field treatment of the topological interaction supports 
this conclusion 

Si. 

Recently, anyons with non-abelian statistics have been the 
subject of great interest because of their potential application 
to topological quantum computation [9]. For these particles, 
there is a degeneracy of the ground state for given positions 
of the particles (so long as they are sufficiently far apart) and 
adiabatically exchanging particles causes the quantum state to 
rotate in this degenerate subspace. Non-abelian anyons are 
hypothesized II lU ll ll 1 121] to exist at the v = | fractional 
quantum Hall state Ill3lll4ill5ll . The vortices of ap x + ip y su- 
perconductor can also have non-abelian statistics, which may 
be observable in Sr 2 Ru04 Ala, Il7l Il8ll . What is the ground 
state of a finite density of non-abelian anyons in zero mag- 
netic field? The general argument about the frustration of ki- 
netic energy applies equally well to non-Abelian anyons. If a 
collection of non-Abelian anyons is bosonic, then it can con- 
dense, potentially lowering the kinetic energy. However, the 



mean-field argument for Abelian anyons does not apply to the 
non-Abelian case. Furthermore, in order for this possibility to 
make sense, the motion of the non-Abelian anyons must break 
the large degeneracy of multi-quasiparticle states. In this pa- 
per, we give arguments which show that, at least in the case 
of non-abelian anyons described by SU (2) K Chern-Simons 
gauge theory, a finite density of non-abelian anyons has a su- 
perfluid ground state. In other words, we expect that doping 
a non-Abelian topological phase in this class has the effect 
of driving a system into a superconducting phase. We note 
that the case of non-Abelian anyons coupled to both SU (2) 
and U(l) Chern-Simons fields has also been shown to be su- 
perconducting when the gauge symmetry is broken down to 

u(i) M- 

Non-Abelian Chern-Simons Anyons. The fractional statis- 
tics of anyons can be implemented in the boundary conditions 
of their wavefuntions, in which case the wavefunctions will 
be multivalued as a function of particle positions. The statis- 
tics can also be implemented by working with single-valued 
wavefunctions (i.e. bosonic or fermionic) with topological in- 
teractions between particles, and the statistics appear in the 
adiabatic transport of particles around each other 
the abelian case, the topological interaction is introduced via 
a fictitious t/(l) gauge field with a Chern-Simons (CS) action 
interacting with the particles The gauge field effectively 
attaches fictitious magnetic flux to each particle, such that the 
wave function accumulates an Aharonov-Bohm phase when- 
ever one particle winds around another. 

Non-abelian fractional statistics can be similarly realized 
by particles interacting via a non-abelian gauge field with a 
CS action: 



S cs = ^J d 3 xe^ x tr (A„d v A x + ^A U A ; 
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where A^ = A a ^T a , T a are the generators of the gauge group, 
and k is an integer, referred to as the level of the theory. We 
limit our discussion to the gauge group SU(2). With each 
particle interacting with this gauge field there is an associ- 
ated representation of the group, its 'isospin'. This internal 
quantum number is totally distinct from the actual spin of the 
particle. The representation carried by a particle or collection 
of particles is physically meaningful, but a particular isospin 
direction can be changed by a gauge transformation. On a 
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closed surface, the total isospin of all of the particles must be 
zero. However, even for given positions of particles in this 
theory, there is a non-trivial ground state degeneracy, corre- 
sponding to different ways in which the isospins can combine 
(or 'fuse') to zero. The allowed isospin states are a truncated 
version of the usual SU (2) spin addition rules, and the trunca- 
tion depends on the level k. Adiabatically moving one particle 
around another rotates the quantum state of the system in this 
space of allowed fusion channels. Hence, the particles have 
non-abelian statistics. Such non-abelian Chern-Simon anyons 
(NACS) appear as excitations of spin models in topological 
phases 0[21, 22, 23]. Therefore it is interesting to ask what 
is the ground state of a system with a finite density of such 
particles that are not static and are able to move around, a sit- 
uation that can arise from doping spin systems in topological 
phases. For the sake of concreteness, we assume the particles 
are in the 'isospin' 1/2 representation. 

Large-K Limit. The Hamiltonian of the CS gauge field 
vanishes, since the Lagrangian in (Q~|) is linear in time deriva- 
tives, but the gauge field itself does not vanish and is deter- 
mined by the matter fields. In second quantized form the full 
Hamiltonian is 



H 



1 

2m 



d 2 r Diip 
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where Di — di — iAf T a , complemented with the constraints: 

\&F% = ..Jh.V: + \^r hc A\A) = -V (3) 
2 J J 2 k 

T a , a = 1, 2, 3 are the generator of SU(2), which are pro- 
portional to the Pauli matrices a a in the spin 1/2 represen- 
tation. The matter field ip(r,t) is a fermionic spin-1/2 two- 
component spinor. 

We consider large n, assuming Af is analytic in 1/n. To 
zeroth order, the A's vanish; and to first order, Gauss' law is 
linear and reduces to the abelian form: 



(4) 



This can be easily solved in the gauge V • A a = 0, similar to 
the abelian case 111 ill . The Hamiltonian to first order in - is: 

H =H + Hi 
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d 2 r^(A a -iV)T a t/j 
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The interaction is weak in the large k limit, and we may use 
the Renormalization Group (RG) scheme of Shankar 112411 to 
analyze it. The RG focuses on the fermionic modes that are 
near the Fermi surface and examines the flow of interaction 



strength in the different channels as the high energy modes 
are decimated. The one-loop flow reveals that only the BCS 
interaction is relevant if it is attractive in any of the angular 
momentum channels. We may write the BCS part of the in- 
teraction in (0 in terms of its angular momentum components 
(suppressing the isospin indices and the momentum sums) : 

2tt i k x k> t t 
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where 9 is the angle between k and k', and A = 
| [kjk' + k' /k) > 1. Depending on the sign of k, the in- 
teraction will be attractive in either the positive angular mo- 
mentum channels or the negative one. Therefore it is relevant 
and there is a BCS pairing instability. This instability justifies 
using a BCS paired wavefunctions as a variational guess for 
the ground state wavefunction of the system, in which fermion 
bilinears (tf)j, a ip^_ k «) develop a non zero expectation value. 
The variational parameters are the isospin structure and angu- 
lar momentum dependence of these bilinear expectation val- 
ues. 

Following the usual BCS approach and assuming the or- 
der parameter has a definite angular momentum dependence, 
Afc = \Ak\e ll ^ k , and either a singlet or a triplet isospin struc- 
ture, we obtain a self consistency gap equation for the order 
parameter: 



a 



l -k'dk' 



-k'dk' 
(6) 



The angular momentum / must be even for spin singlet pair- 
ing and odd for spin triplet, and the interaction is effectively 
stronger in the singlet channel C s — 3Ct- We find that the 
lowest energy paired state is one with isospin singlet pairing 
and I = 2 pair angular momentum. Since an isospin singlet 
does not couple to the SU(2) gauge field, such a pair is a bo- 
son. If the anyons are charged, this is a superconductor with 
d + id pairing. 

One might worry that the BCS approach is not valid be- 
cause the topological interaction between quasiparticles is 
long-ranged. However, when two SU(2) K anyons pair to 
form a boson, the resulting bosons do not have any remaining 
long-ranged interactions. By pairing in singlets, the anyons 
lose their non-abelian interaction, and the ground state can be 
thought of as condensation of pairs, similar to condensation 
of pairs of fermions in superconductors. Therefore, the long- 
ranged part of the topological interaction is not important in 
this finite-density situation. At distances longer than the su- 
perconducting coherence length, the topological interactions 
do not affect the important degrees of freedom. 

k = 2 Case. We now argue that a phase with excita- 
tions satisfying the SU (2)2 NACS statistics can be reached by 
quantum disordering a superfluid, which implies the inverse: 
k = 2 NACS anyons can condense into a superfluid. 
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To disorder a superfiuid, it is useful to invoke a duality 
transformation and write the effective action in terms of vortex 



variables [25, 26, 31]: 



S„ = / drd 



(|(5 / ,-io„V| a + WI 2 )+ 



Q Z7T 



(7) 



where ip v is the vortex annihilation operator. The vortices 
interact logarithmically, as expected, via the gauge field a M , 
which is related to the super-current by: = ^e^ v \d v a\ 
and A em is the external electromagnetic potential that cou- 
ples to the charge current. The superfiuid and insulator phases 
can be recovered from this effective action: When the vor- 
tices are gapped, ((p v ) = 0, the vortex field and the gauge 
field can be integrated out and the resulting effective ac- 
tion contains a mass term for the external electromagnetic po- 
tential, i.e. the system is a superconductor. When vortices 
condense, (ip v ) ^ 0, the resulting effective action for A em 
has a Maxwell term to leading order, so the system is insu- 
lating. Condensation of any powers of the vortex field would 
also result in an insulator. 

Recently, it has been argued that it is justified in certain lim- 
its to treat vortices as weakly interacting fermions rather than 
strongly interacting bosons 132113311 . Fermionizing is achieved 
by attaching one Abelian Chern-Simons flux quantum per par- 
ticle. This procedure makes vortices fermionic at the expense 
of adding a gauge field with a CS action: 

A> = -07 M (du - «<V - ia^ip + V{^) 

+ TT^^^^ax) 2 + —e^xa^ax + -^-A^ju 
Zg 47T Ztt 

(8) 



where ip is the usual two-component spinor of Dirac fermions 
in 2 + 1 dimensions and the 7 M can be chosen to be Pauli ma- 
trices. In the presence of a gauge field with a Maxwell action, 
it is possible to integrate out the CS field and this will gener- 
ate only terms that are higher order than the Maxwell term for 
the remaining gauge field; naively, they are less relevant. This 



is done by shifting a p 



a p to decouple a M from the 



fermions. Then can be integrated out to obtain: 

C v = ftf(d u - iajtl) + m(^0 + G(^) 2 

+ ^(e^xd.ax) 2 + + 0(d 3 a 2 ) (9) 

where we have expanded V(ipip). We will be interested in 
a phase in which spontaneous symmetry breaking leads to a 
mass for the gauge field a^, and therefore higher-order terms 
in a M can be treated perturbatively. We ignore them here. 

In terms of fermionic vortices, an insulating phase is 
reached when pairs of vortices condense. Pairing of fermions 
with a relativistic spectrum has been explored in two 02911 and 
three BOIl spatial dimensions. We can decouple the quartic 



interaction term in (0 via a Hubbard-Stratonovich transfor- 
mation and look for a self consistent pairing Ay : 



Sf = J drd 2 xipj ll (d u — ia^ij} + rmpip+ 
+ %ft (x)iAij {x)ift (x) + h.c. 



(10) 



It was found M29I1 that for G > no such pairing exists, while 
for G < the only self consistent pairing is A = j 2 A(x), 
where A(x) is a scalar. When such pairing happens, the sys- 
tem is in the electrically insulating phase. 

In the presence of pairing, the mean field Hamiltonian is 
quadratic and is diagonalized by fermionic operators, ajk, that 



are linear combinations of ipi an d V4 1 3 



(11) 



the i and j indices are component indices, and the index k 
enumerates the state, for example it can be momentum in a ho- 
mogeneous system. The functions Uijk (x) and Vijk (x) solve 
the Bogoliubov deGennes (BdG) equations: 
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V «2(*0 J 
(12) 

with ho = 7°(— iV-7+m) the Dirac Hamiltonian. Equations 
of similar structure have been discussed in ref. [34] in the con- 
text of graphene, with the important difference that the order 
parameter there corresponds to a mass rather than pairing of 
the fermions. 

The pairing A(x) is much like the familiar pairing of non- 
relativistic fermions. Most importantly, it gives the dual gauge 
field a mass, and the system is a "dual superconductor", i.e. 
an electric insulator. The effective action for A{x) is simply 
a Landau-Ginzburg free energy functional of a three dimen- 
sional superconductor, and we will draw on our knowledge 
of the conventional theory of superconductivity to learn about 
the phase described above. 

The boson charge density maps under duality to the dual 
magnetic flux j — V x a. so changing the charge density 
is equivalent to introducing dual magnetic flux to the system. 
In the insulating phase where there is vortex pairing, the dual 
magnetic flux can form a lattice of defects in the vortex pair 
order parameter A(x), each carrying half a unit of dual flux, 
i.e. an Abrikosov lattice. This configuration corresponds to 
a crystal of localized charges, each carrying half the charge 
of the original bosons which formed the superconductor. If 
the bosons are charge 2e Cooper pairs, then the paired-vortex 
insulating phase has localized excitations of charge e. We 
will now show that these excitations also obey SU (2)2 NACS 
statistics. 

A prerequisite for having non-abelian statistics is having 
a degeneracy for a given position of the excitations. In the 
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paired-vortex insulator described above, the degeneracy is due 
to zero modes localized inside the defects in the pairing field 
A(x). For simplicity, a defect can be modeled as a region 
where the magnitude of the pairing field A(x) drops abruptly 
to zero and its phase changes by 2tt in going around the defect, 
A(r, 0) = 9(r - r )A oo e ±ie with 0(s) the step function, 
and ro is on the order of the coherence length of the conden- 
sate, £„. Solving the BdG equations, eq. (fT2b . with this A(r, 0) 
gives the spectrum of excitations in the presence of a defect. 
We are interested in solutions at zero energy, and if we also 
restrict ourselves to solutions satisfying u. L (x) = v*(x), then 
the BdG equations reduce to only two linearly independent 
equations: 

mm + e-' w (d r - -d e )u 2 + A(r)e ±ie u* 2 = 
r 

-e l9 (d r + -d e )m ~ mu 2 - A(r)e ±l %* = (13) 
r 

we can assume m > and Aoo > without loss of general- 
ity. These equations can be solved exactly, but it is enough to 
examine their behavior for large r. The asymptotic behavior 
of m an d u 2 is: 

ui, 2 (r -► oo) ~ Ae-^ +m '> r + B e - {A ™- m)r (14) 

For Aoo < rn, there is no well behaved localized solution, 
but for Aoo > m there is one solution which decays exponen- 
tially away from the defect. The fermionic operator = 7^ 
corresponding to this zero-energy solution localized at the de- 
fect has the property 7^ = 7, i.e. it is a Majorana fermion, 
as can be seen from eq.dTTb with Ui(x) = v* (x). A two-state 
Hilbert space cannot be built from a single Majorana fermion, 
but in the presense of two defects in A(x), which correspond 
to two localized charges that are well separated, there is a Ma- 
jorana mode in each of the defects, which allows us to define 
a Dirac fermion \P = ^(71 + ^72) satisfying the usual an- 
ticommutation rules, ^>'} = 1, The occupation number 
of the fermionic mode W defines two degenerate zero energy 
states, |0) and W^|0). This degeneracy does not correspond 
to a local quantum number since the Dirac fermion is shared 
between the defects, regardless of their separation. For 2N 
defects in the the pairing field A(x) there are 2^ degenerate 
quantum states at zero energy. The non-abelian statistics is 
evident when exchanging positions of defects. When defects 
i and j are exchanged, the Majorana zero modes are trans- 
formed by 7j — * 7j and — ► —7,, which can be translated 
into a rotation in the 2 N dimensional degenerate subspace of 
states JT^B. This type of non-abelian statistics corresponds 
to SU{2) 2 . 

Discussion. In this paper, we have studied the ground 
state of a finite density of dynamical non-abelian SU(2) K 
Chern-Simons anyons. In the limit k — » 00 we showed that 
the topological interaction between anyons induces a pairing 
instability, and the preferred pairing has a d + id spatial sym- 
metry and j — isospin structure. For the special value 
k = 2, we showed that a superfluid of pairs can be disordered 



to get an insulator with SU(2)2 anyonic excitations, which im- 
plies the inverse path is also possible, and SU(2)2 anyons can 
pair to form a superfluid. Since pairs of NACS anyons in the 
j = channel are bosonic with respect to each other for all 
values of K, we suggest that the ground state for a finite den- 
sity of NACS anyons is a superfluid of bosonic pairs for all 
K. 

The most obvious implication is that if a frustrated mag- 
net were found in a non-Abelian topological phase, then dop- 
ing the system would drive it superconducting. Conversely, 
non-Abelian anyon superconductivity may be an explanation 
for the superconductivity of some anomalous materials. For 
the route to superconductivity which we sketched above for 
k = 2, vortex pairs condense en route to the non-Abelian 
insulating state. Consequently, individual vortices do not con- 
dense. Therefore, flux-trapping experiments 1 36] should find 
a positive result for such a superconductor. 
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